Universal, additive effect of temperature on the rheology of amorphous sohds 



O 



i> 



Joyjit Chattoraj'^^', Christiane Caroli*^^', and Anacl Lcmaitre'-^^ 

'^' Universite Pans Est - Lahoratoire Navier, ENPC-Paris, LCPC, 

CNRS UMR 8205 2 allee Kepler, 77420 Champs-sur-Marne, France and 

'■^' INSP, Universite Pierre et Marie Curie-Paris 6, CNRS, 

UMR 7588, 140 rue de Lourmel, 75015 Paris, France 

(Dated: May 10, 2010) 

Extensive measurements of macroscopic stress in a 2D LJ glass, over a broad range of temperatures 
(r) and strain rates (7), demonstrate a very significant decrease of the flowing stress with T, even 
much below the glass transition. A detailed analysis of the interplay between loading, thermal 
activation, and mechanical noise leads us to propose that over a broad (7, T) region, the effect 
of temperature amounts to a mere lowering of the strains at which plastic events occur, while 
the athermal avalanche dynamics remains essentially unperturbed. Temperature is then shown to 
correct the athermal stress by a (negative) additive contribution which presents a universal form, 
thus bringing support to and extending an expression proposed by Johnson and Samwer [1]. Our 
prediction is shown to match strikingly well numerical data up to the vicinity of Tg. 
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It is now well accepted that macroscopic plastic defor- 
mation of amorphous solids is the net result of an ac- 
cumulation of local rearrangements ( "shear transforma- 
tions" , or simply "flips" ) involving small clusters of atoms 
or particles [2]. However, what controls the occurrence of 
these flips, ie when and where they take place, remains a 
controversial question, which needs to be answered prior 
to the formulation of any theory of plasticity. 

Recent numerical works bring fragmentary informa- 
tion about these mechanisms. In athermal simulations- 
whether quasi-static or at finite strain-rate-it is possible 
to identify, at any time, in the steady-state sheared sys- 
tem a set of soft zones, which retain their identity for 
sizeable stretches of strain [3]. As a zone is loaded by 
the external drive, it gradually softens until it reaches 
its instability threshold and flips [4]. This event creates 
in the surrounding medium a long-ranged elastic fleld, 
with quadrupolar symmetry [5], which shifts the strain 
of other zones, hence may induce secondary events. This 
mechanism gives rise to avalanches of flips [6] with a 
strain-rate dependent average size [7] . 

These results have led us to propose a picture in which, 
at T = 0, the dynamics of soft zones in the vicinity of 
their instability thresholds plays a critical role in control- 
ling dissipative events, while structural disorder, which 
permits the existence of soft zones in the flrst place, is 
otherwise accessory to the unfolding of plastic events. 
But we must then ask what is the effect of a flnite tem- 
perature on plastic activity [8]. Indeed, recent works 
concerned with stress and elastic fluctuations, as well as 
energy barrier distributions, have shown that strained 
systems under steady flow present a sizeable fraction of 
low-energy barriers, most of which are not involved in 
the athermal response [9]. It is then possible that, even 
a small temperature may activate jumps over these many 
other barriers, hence disrupting the avalanche processes 
which control the athermal response. 

Here, we perform extensive measurements of the 



macroscopic stress in a 2D Lennard-Jones glass, over a 
broad range of temperatures (T) and strain rates (7). 
We find that finite temperatures induce a very signifi- 
cant decrease of the fiowing stress. We show that these 
data can be interpreted within the framework provided 
by a detailed analysis of the interplay between loading, 
thermal activation, and mechanical noise. This analysis 
predicts that over a broad region of the (7, T) parameter 
space, the avalanche dynamics should remain basically 
unperturbed, the effect of temperature amounting to a 
mere lowering of the strains at which plastic events oc- 
cur, leading to an expression for stress of the form: 

a(7,r) = 'To(7) + $(7,T) 

where (i) cro(7) exhibits precisely the dependence found 
in athermal systems and (ii) the (negative) shift func- 
tion $ ex -(-rin(C7/T5/6))2/3 is analogous to that 
proposed by Johnson and Samwer [1] to account for a 
number of experimental results on metallic glasses. This 
expression matches strikingly well our numerical data up 
to the vicinity of the glass transition temperature Tg. 
We thus conclude to the robustness of avalanches, which 
also show up clearly in the same temperature range in 
the structure of strain maps. 

We use the same 2D binary LJ mixture as in Refs. [3, 
7]: it is composed of large (L) and small (S) parti- 
cles with radii Rl ~ 0.5, Rs = 0.3 and equal masses 
m ~ 1 (in standard LJ units), in a number ratio of 
Nl/Ns = (1 -I- •\/5)/4 to ensure that no crystallization 
occurs, at packing fraction 'k{NlR\ + NsR%)/L'^ = 0.9. 
With these parameters the shear modulus is ~ 20 [3], 
the shear wave speed Cs — 3.4, and the glass transition 
temperature (identified as the value at which Tq = 10^, 
and computed following [10]) is ^ 0.27. Finite temper- 
ature MD simulations are performed on square L x L 
systems, using Lees-Edwards boundary conditions, with 
a standard velocity rescaling protocol [11] and time steps 
dt < 0.01. All the data presented here are obtained af- 



ter 100% of strain, which ensures that the system is in 
steady state. Large sets and long strain intervals are used 
for statistical accuracy (e.g. for L = 40, 25 samples are 
strained up to 1300%). 

We present on Fig. l-(a) steady state stress a vs strain 
rate 7 for different system sizes L = 10, 20, 40, 80, 160, 
two temperatures T = 0.025 and T = 0.2, and 7 ranging 
from 4 X 10~^ to 10~^. For each temperature we find, as 
in athermal systems, quick convergence with increasing 
L, of (1(7) towards a master curve. In the range of 7 
studied, saturation is already reached for L = 40, which 
allows us to focus in the following on stress data obtained 
for this system size. 
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FIG. 1. (a): The macroscopic stress cr as a function of strain 
rate 7 for systems of sizes L — 10, 20, 40, 80, and 160, at 
temperatures T — 0.025 and T = 0.2. (b): a versus 7 for 
L = 40, compared with fits of the form a = ^0 + ^1 VT" (solid 
lines). 



To qualify the effect of temperature, we then attempt 
(see Fig. l-(b)) to fit these curves with a = Aq + Ai ^^, 
a form which was shown to match very well the results of 
athermal simulations. For our lowest temperature T ~ 
0.025, the fit is quite satisfactory, although not as good 
as was found with athermal data. This suggests that at 
T finite but <C Tg the dynamics is closely similar to that 
of an athermal system. At higher temperature, however, 
not only does this fit become definitely poorer but, more 
importantly, the overall stress level drops significantly: 
temperature thus has a pronounced effect on dissipation. 
Does this mean that it completely modifies the avalanche- 
dominated dynamics at work in the T — > limit? 

We know that plastic events correspond to local shear 
transformations (zone flips), which are preceded by the 
gradual vanishing of an energy barrier under loading. 
At T = 0, each flip is triggered right at the threshold 
strain jc where instability is reached. At finite T, we 
must expect that flips occur in anticipation due to ther- 
mal activation, which becomes gradually more efficient 
when a zone approaches its 7c [12]. More precisely, con- 
sider a single zone, initially lying at a strain 70 < 7c, 
which is loaded at finite strain rate 7 in the presence 
of thermal noise. At an external strain 7 € [7o,7c[, 
the probability ^(7) that it has not yet flipped obeys 



P(7 + d7) = P(7)(l-aii?(7)),or: 

— ==--i?(7)P(7;7o) 
with R{"f) the rate of activated jumps. The solution is: 



-P(7;7o) =exp 
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At low temperature, activation is efficient only close to 
the saddle node bifurcation occurring at 7c. There, the 
energy barrier and the "attempt frequency" present the 
universal scaling forms: AE = B {'fc — 7)'^^^ and u) = 
^{ic - if'^ [4] [13]. Provided that AE/T > 1, we can 
use the standard Kramers expression for the activation 
rate: i?(7) = w exp(-A£;/r), so that: 



-P(7;7o) = exp 




Q{h) - ^(-570) 



where J7 = 7c - 7 and Q{5-i) = F (|; f 5^1^/'^), with 
F the upper incomplete gamma function. Our usage of 
Kramers' expression implies that e = T/B5^'^''^ <ti 1, 
so that we can use the asymptotic expression F(s; x) ^ 
x"^^ e^^, whence Q{Sj) ^ e^^^ exp(— 1/e). 
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FIG. 2. The parameters used come from from the fit of our 
stress data (see text). For T = 0.1: (a) the function P vs 
57 for 70 — —00 (solid line) and 70 ^ 7* (dashed); (b) 5^ '^ 
vs 570 (sohd line) with asymptotes (thin hue), (c): 7''™' and 
•y^''' vs T (using a/d — 1/3, Aeo = 5%); shaded area: the 
high-T region where data start to depart from fit. 

Plots of P(7;7o) are presented on Fig. 2-(a). Let us 
first consider the formal limit 70 — >■ —00 {Q{5^o) — > 0), 
which corresponds to the case when a zone is initially 
very far below threshold. Because P is an exponential of 
an exponential it presents a very sharp transition from 
P ~ 1 to P ^ around a strain 7* such that: 



li(l)'%<-^-)- 



which is next solved at leading order in e, yielding: 
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The width of this transition 
|9P/97(7*;-oo)|-i = (2e/3)e*,57* 
0{e*) relative to Sj* itself. 

As 7o increases from — oo, the curves P(7;7o) remain 
nearly identical to P(7; — oo) up to the immediate vicin- 
ity of 7"^. For 7o > 7"^, P(7;7o) no longer presents a 
plateau at low values of 7, but drops already sharply at 
70, with slope \dP/d-f{-fo] 70)] = i -R(7o), which is an in- 
creasing function of 70. Therefore, the solution 7^'P of 
P(7^'P; 7o) — 1/e can always be interpreted as the typical 
strain at which a zone starting from 70 flips. 

The curve 7*''p(7o) (Fig. 2-(b)) exhibits a sharp tran- 
sition between two limiting behaviors: (i) when injected 
at 70 ^ 7*, a zone flips at 7^^'^ ~ 7*; (ii) when injected 



at 70 
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it flips almost immediately: 7^'P 
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The width of the transition can be estimated as be- 
ing 0{5^* €*) ^ 57*, which compares with the width 
of the drop of P(7; — 00) around 7*. We can therefore 
conclude that the competition between thermal activa- 
tion and drive defines (slightly fuzzy) apparent thresh- 
olds which are shifted by — 57*(7, T) from the mechanical 
yield points. 

In our previous studies of athermal systems, we found 
that the mechanical noise-that is the stress noise gener- 
ated by the flips themselves-played a key role by inducing 
correlations between flip events, leading to the emergence 
of avalanche behavior. At finite T, a zone embedded in a 
sheared system is thus experiencing both thermal noise 
and random strain shifts (ie barrier height fluctuations) 
originating from prior events. In [7], we proposed to sepa- 
rate mechanical noise into: (i) a low-frequency part, gen- 
erated by nearby events (within a sphere of radius €); (ii) 
a background noise coming from the rest of the system. 
Requiring that the near-field signals be non-overlapping 
and stand out of the background noise defines a single 
length £(7) which we identified as the avalanche size. 

An exact treatment of the jump dynamics in the pres- 
ence of activation, loading, and barrier fluctuations is 
for the moment out of reach. However we note that 
the values 7* can still be interpreted as effective thresh- 
olds, provided that the unfolding of the activated jump 
of a zone near its 7^'P is not perturbed by the ambi- 
ent noise (which includes nearby and far- field signals). 
For this purpose, we recall that the strain field due to 
a shear transformation has the form /esh. = ^Aeq [5], 
where a is a zone size and Aeq an elementary plastic 
strain. In a system of size L, the flip rate is 7?. == 
7L^/(a^Aeo), and the strain fluctuation due to the noise 
incoming from the whole system during a time r ver- 
ifies A72(r) = l^jj^fl^d^ = 7t^, with d a 
typical inter-zone distance [7]. In order for the activa- 
tion process to be negligibly perturbed, we must ensure 
that the strain fluctuation accumulated during the acti- 
vation time l/i?(7^'P) < 1/^(7"*), remains much smaller 
than the thermal strain shift 5^* itself. This is guar- 



anteed as soon as: j a^Aeo/{R{j*)cr) <ti (Sj*) . As 
-R(7*) = It 7 V^T*' tliis ^Iso writes: 
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When this condition is satisfled, the process of flip 
activation disentangles from the response to incoming 
mechanical noise signals. The basic elements of a phe- 
nomenology of avalanche dynamics, namely the advcc- 
tion of zones towards (shifted) thresholds and the pres- 
ence of mechanical noise signals, are preserved. In par- 
ticular, the separation between nearby, correlated, signal 
and background noise proceeds exactly along the same 
lines as in the T = limit [7], thus defining the same 
value for the avalanche size ^(7). As flip events occur at 
shifted thresholds, the macroscopic stress should thus be 
of the form: 



cr(7,r) = 0-0(7) ~ t^S-f*{^,T) 



(5) 



where /i is the shear modulus. From athermal simu- 
lations [7], we know that uo is of the form: cto(7) = 
Aq + Ai y/^. The average (S7* accounts for the fact that 
the values of B and v appearing in the above calculation 
are distributed due to structural disorder. Assuming that 
In B and In v are well-centered, we can finaly write: 



a(7) ^Aa+Ai^-A2 P'/' [in (a^ T''/^/^ 



2/3 



(6) 
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with A2 = ^P~2/3 and A3 = 

To test this prediction, we now fit the stress data shown 
on Fig. 1 using this four parameter expression. As seen 
on Fig. 3, the fit is remarkable over quite a broad range 
of parameters, T varying from nearly up to ~ 0.2 (to be 
compared with Tg = 0.27), and 7 ranging over more than 
two decades. The values of the parameters are strongly 
constrained by the fit: Aq and Ai by the low-P data; A2 
by the low-7 data. It clearly confirms the form of cto (7) as 
well as the functional form of the correction term, except 
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FIG. 3. Macroscopic stress a (filled symbols) as a function 
of strain rate 7 (left) and T (right) compared with the fit 
obtained using Eq. (6), with parameters; Ao — 0.66, Ai = 
2.09, A2 = 0.27, and A3 = 0.22. 



for the 5/6 exponent which is weakly discriminated. The 
T^' ■^ dependence, which is very visible at low 7, is a clear 
signature that the saddle-node bifurcation controls the 
behavior of energy barriers heights over the whole range 
of relevants strains. 
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FIG. 4. The non-affine strain field accumulated over strain 
intervals A7 — 1%, 5%, 10% (left to right), when sheared at 
7 = 10"'', and for T = 0.025, 0.2, 0.3 (bottom to top). System 
size is 160 x 160. 



From the fits we extract typical values v ^ 50 and 
B ~ 650, from which we can examine a posteriori the va- 
lidity of our expression for stress. The applicability of the 
Kramers expression (1/e* = ln((2t//37) {T/Bf/'^) > 1) 
and condition (4) both set upper limits 7^''' and 7^'™- 
(resp.), which are increasing functions of T (see Fig. 2- 
(c)). By iteration, we have used for our fit only the 
points which satisfy these conditions. This excludes a 
few data points at the lowest temperatures and highest 
strain rates. 

The conditions 7 < i^'^T), 7 < ^^''^■(T) and T < 0.2 
define a large region of parameter space within which we 
claim that the effect of thermal fluctuations reduces to a 
universal additive contribution to stress of the form given 
in Eq. (5). In the vicinity of Tg, the gradual departure 
of the measured stress away from this expression should 
be assigned to the increasing contribution of thermally 
activated events of a different nature than the zone flips 
that we have considered here. We believe, these might be 
related to the finite, low-lying barriers shown by Rodney 
and Schuh [9] to be present in sheared systems. 

The present work provides a firm support to the notion 
of effective threshold and generalizes significantly the ex- 
pression shown by Johnson and Samwcr to fit a large 



body of experimental data on metallic glasses [1]. In- 
deed, our derivation clarifies why an additive correction 
to stress, with the universal form set by the saddle-node 
bifurcation, holds despite the presence of structural dis- 
order. Moreover, our data show that this correction adds 
up, as predicted by our model, to a stress 0-0(7) which is 
precisely that expected in an athermal system. We thus 
conclude to the robustness of avalanche dynamics up to 
the vicinity of Tg . 

In order to gain further insight into this question, we 
now show how the non-affine strain field accumulates, 
at different temperatures, as the system is macroscopi- 
cally sheared. Strain maps are displayed on Fig. 4 for in- 
creasing values of the external strain A7 = 1%, 5%, 10% 
(left to right), and for T = 0.025,0.2,0.3 (bottom to 
top). Strikingly enough, at all values of T, we clearly 
see comparable correlated structures. The strain pat- 
terns at T = 0.025 and 0.2 are similar and exhibit, as 
those seen previously in athermal simulations, a clear di- 
rectionality which can only be due to Eshclby-likc clastic 
interactions [7]. It is only at T = 0.3, i.e. in the super- 
cooled regime, that a qualitative change can be observed, 
the structures becoming slightly blurred and shorter- 
ranged. We thus see a continuity of behavior across 
the glass transition which connects the low-temperature 
avalanche behavior with the cooperative dynamics ob- 
served by Tanaka in the supercooled regime [14]. More 
quantitative information on this question might come 
from the detailed analysis of the (7, T)-dependence of 
the diffusion coefficient, a study which is presently under 
way. 
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